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ABSTRACT 

We study the star forming regions in the spiral galaxy NGC 432 1 (M100). We take advantage 
of the spatial resolution (2.5 arcsec FWHM) of the Swift/UVOT camera and the availability 
of three UV passbands in the region 1600 < A < 3000A , in combination with optical and 
IR imaging from SDSS, KPNO/Ha and Spitzer/IRAC, to obtain a catalogue of 787 star form- 
ing regions out to three disc scale lengths. We use a large volume of star formation histories, 
combined with stellar population synthesis, to determine the properties of the young stellar 
component and its relationship with the spiral arms. The Ha luminosities of the sources have 
a strong decreasing radial trend, suggesting more massive star forming regions in the cen- 
tral part of the galaxy. When segregated with respect to NUV-optical colour, blue sources 
have a significant excess of flux in the IR at 8/xm, revealing the contribution from PAHs, al- 
though the overall reddening of these sources stays below E(B-V)=0.2 mag. The distribution 
of distances to the spiral arms is compared for subsamples selected according to Ha lumi- 
nosity, NUV-optical colour, or ages derived from a population synthesis model. An offset is 
expected between these subsamples as a function of radius if the pattern speed of the spiral 
arm were constant - as predicted by classic density wave theory. No significant offsets are 
found, favouring instead a mechanism where the pattern speed has a radial dependence. 

Key words: galaxies: spiral - galaxies: star formation - galaxies: general - galaxies: individ- 
ual (NGC 4321) - 
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1 1 INTRODUCTION 



X 



Ninet y years after the Great De bate on the the nature of "spiral neb- 
ulae" dShaplev&CurtIslll92ll) , our knowledge about disc galax- 
ies has grown considerably. However, the processes responsible for 
their spiral appearance are still under scrutiny. Various mechanisms 
have been proposed, the widely accepted one being the triggering 
of a long-lived den sity wave that mo ves as a solid body, with con- 
stant pattern speed dlin & Shul 19641) . One can define the corotation 
radius as the point at which gas and stars move with the same speed 
as the spiral pattern. Since the rotation velocity of gas and stars is 
roughly independent of radius, the regions inside (outside) the coro- 
tation radius will move faster (slower) than the spiral arm pattern. 
The gas component piles up in th e spiral arms, e xperiencing a shock 
that gives rise to star formation dRobertsll 19691) . One would expect 
in this scenario that the youngest stars bora from the molecular 
clouds in the spiral arms would be found slightly ahead of the arm 
traced by the molecular gas, if located within the corotation radius, 
and behind the arm outside of corotation. Therefore, if we observe 
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the tracers of the different stages of star formation, such as HI, CO 
and Ha, we expect to find an offset among them as a function of 
the radius. There are a nu mber of observational efforts to detect 
such offset in the past (e.g. IVogel et al . 1988; Garcfa-Burilloetal 



1993; 



. Randl i 

Egusa et aU2009l : lFovle et alj|201 lj). The majority of these studies 



1995), and more recen tly (e.g 



Tamburro et al. 200 



measure the angular offsets between CO as a tracer of star forming 
molecular clouds and Ha and/or 24/im as a trace r of young stars. A 
time difference around a few Myr is found (e.g. lEgusa et al]|2009l . 
and references therein) between CO and Ha. We note that HI is 
not necessarily a direct tracer of molecular cloud formation, as it 
often results from molecular gas dissociation caused by star forma- 
tion activity dAUen et al.ll98 6: Young & Scovillell99ll ). A stronger 
constraint can be imposed if we target instead tracers that cover a 
wider range of ages for the distribution of recently formed sources. 
The use of UV light opens up this possibility, since s tellar popu- 
lation s can be tracked over the first few hundred Myr dFovle et"al] 
1201 lh . 

An alternative mechanism for the triggering of spiral arms 
involves self-gravitat ional instabilities (Goldreich & Lvnden-Belll 
1 19651 : iToomrel I 198 lh whose spiral features have much shorter 
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Table 1. Log of SwiMJVOT observations of NGC 4321 used in this paper. 



Exposure Time (s) 



ObsID 


Date 


UVW2 
1928A 


UVM2 
2246A 


UVWl 
2600A 


00035227001 


2005/1 1/06 


532.0 





80.6 


00035227002 


2005/11/08 


386.6 


290.1 


192.7 


00035227003 


2005/11/13 


6832.3 


5107.6 


3414.5 


00030365002 


2006/02/08 


673.6 


505.6 


332.2 


00030365005 


2006/02/10 


1023.6 





343.1 


00030365006 


2006/02/12 


1092.1 


390.6 





00030365007 


2006/02/14 


939.6 


704.6 


462.1 


00030365008 


2006/02/17 


998.0 


745.0 


486.2 


00030365009 


2006/02/18 


17.8 








00030365011 


2006/02/20 


595.4 


298.5 


98.4 


00030365012 


2006/02/21 


325.6 


244.6 


160.9 


00030365013 


2006/02/24 


543.2 


405.2 


266.8 


00030365014 


2006/02/26 


519.2 


387.2 


253.0 


00030365015 


2006/02/28 


769.2 


578.2 


377.9 


00030365016 


2006/03/02 


443.4 


332.4 


159.9 


00030365017 


2006/03/05 


701.3 


505.6 


331.2 


00030365018 


2006/03/12 


1032.4 


773.4 


508.2 


00030365019 


2006/03/18 


1103.4 


827.4 


542.7 


00030365021 


2006/03/27 


1045.7 





468.1 


00030365023 


2006/03/29 


1780.2 





1368.0 




TOTAL (ks) 
Effective TOTAL (ks) 


21.35 
21.28 


12.10 
11.95 


9.85 
9.74 



lifetimes. Numerical simulations have not been capable of ex- 
plaining the onset of a stable density wave with constant pattern 
speed, whereas short-lived spiral arms appear natura l ly (see e.g. 
Sellwoodl201 it iFuiii et al.ll201 it lOuillen et alj|201ll ; IWada et al.1 
201l tlGrandetal.l2012l) . Other processes, like tidal interactions or 
star-formation induced structures have been proposed as well, leav- 
ing a differe nt imprint on the prope rties of the underlying stellar 
populations jDobbs & Pringlell2010l) . 

In this pap er we make use of the Ultraviolet/Optical Tele- 
scope (UVOT, Roming et al. l2005t) on board the Swift space- 
craft dGehrels et al.ll2004h . to study in detail the star forming re- 
gions of a nearby grand-design spiral. NGC 4321 (M100) is a 
face-on late-type spiral gala xy (SABbc) towards the Virgo cluster 
dde Vaucouleurs et al]|l99ll) . Even though Swift is an observatory 
built for the study of transient phenomena, most notably Gamma 
Ray Bursts, the resolution of UVOT makes it an optimal instrument 
to study star formation in nearby galaxi es. The presence of a recen t 
supernova in NGC 4321 (SN2006X, ISuzuki & Migliardil 120061) 
triggered a number of follow-up observations with Swift/UVCfT, 
allowing us to put together one of the deepest UV images of a 
nearby galaxy at a spatial resolution of 2.5 arcsec (FWHM) over 
a ~ 100 arcmin 2 field of view. In this paper we adopt the dis- 
tance from the Hubb le Space Telescope Key Project to measure 
Hubble's constant dFreedman et alJl200ll) . with a distance modu- 
lus n = 30.91 ± 0.07 mag, implying a distance of 15.2 Mpc, and 
an angular scale of 73.7 pc/arcsec. The scale length of the disc of 
NGC 4321 is h=1.39 arcmin or 6.15 kpc measured in the V band 
jKodairaetai]|l986h . We use this value as a characteristic disc 
length throughout this paper. The radial positions of the sources 
are deprojected by using an inclination of 38°, with the re ceding 
side located at a position angle of 151 dChemin et ai]|2006l) . 



2 THE UV IMAGES 

We extracted the UV images of NGC 4321 from the archive of 
the Swiff/UVOT missiorQ at HEASARC. The log of the obser- 
vations can be found in Table Q] UVOT is one of the three in- 
struments on board the Swift observatory. It couples a 30 cm 
Ritchey-Chretien telescope to a photon-counting, microchannel- 
plate-int ensified CCD camer a. Details of the instrument can be 
found in iRoming et ai] d2005h. whereas the calibration is described 
in lPoole et al.1 d2008) and B reeveld et alj d2010L l201 ]]). We use the 
three available UV filters: UVW2, UVM2 and UVWl, that span the 
wavelength range: 1600-3400 A . These three filters roughl y map 
the sp ectral window of the GALEX NUV passband (Mar tin et al.l 
l2005h . with UVM2 centered at the position of the 2 175 A UV-bump 
in the Milky Way extinction curve dConrov et alfeoiol) . Fig. 10 in 
IPoole et all d2008l) shows the in-orbit effective area of the UVOT 
filters. The spatial resolution of the instrument, along with a Point 
Spread Function of 2.4, 2.5 and 2. 9 arcsec (FWHM) in UVW2, 
UVM2 and UVWl, respectively dBreeveld et al] 1201 Oh . makes 
UVOT images significantly better resolved than NUV images wit h 
GALEX, which has a FWHM of ~6 arcsec dMorrissev et al]2005l) . 
We note that at the distance of Ml 00, the UVOT resolution limit 
maps a physical distance around 200 pc, which is the maximum 
possible value to be able to neglect deble nding issues in HII re- 
gions, given their characteristic sizes dKennicuttetal.lll989h . 

We retrieved the science-grade exposures from the archive. 
Some of these images did not have an aspect correction. This fail- 
ure is often caused by the lack of nearby stars in the USNO cata- 
logue, that are used as reference in an automatic way to determine 
this correction. In order to have as many images as possible, we 
wrote our own software to aspect correct the images from a man- 
ual selection of foreground stars, and to combine all images with 
a flux-conserving drizzling algorit hm at the same pixel si ze (i.e. 
0.5 arcsec), based on the method of iFruchter & Ho ok (2002). Some 
of the images, with exposure times typically below 30 seconds, 
could not be used because they could not be properly registered. 
Hence, we define as effective integration times those from the fi- 
nal exposure map within the footprint of NGC 4321. The exposure 
map in this region is homogeneous. The total effective exposure 
times of the UV data are 2 1 .28ks in UVW2, 1 1 .95ks in UVM2 and 
9.74ks in UVWl. The top panels of Fig.Q]shows greyscale images 
of the UVOT combined exposures. We illustrate in Fig.|2]the differ- 
ence between our UVOT/UVW2 deep exposure and the available 
archival image from GALEXfNUV (exposure time 1183 s). Both 
frames show the central 8x8 kpc 2 region. 

The source dete ction was performed with SExtractor 
dBertin & Arnoutsll 19961) on the UVW2 image, which is the deep- 
est one of the three, and corresponds to the shortest wavelength 
- optimally mapping ongoing star forming regions. The SExtrac- 
tor code has been optimized for the detection of sources over a 
homogeneous background (e.g. galaxies in a survey or stars in a 
globular cluster). However, our case is slightly different, as the star 
forming sources are located on top of a diffuse background that 
has to be removed before detecting the sources. To contend with 
this, we define our own background on the original UVW2 im- 
ag e, following a p r ocedu re inspired by the methodology set out 
in iHoversten et al.l d201ll) . A number of images are obtained by 
performing a circle median of the observed image over a range of 
radii. In order to determine the background, we need to subtract 
the information over smaller scales, ranging from individual HII 
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The Swift/UVOT view ofNGC 4321 3 




IRAC Sum 



Figure 1. Views of NGC4321 (M100) in different passbands over a 8 X 8 arcmin 2 field of view. North is up, East to the left. The top panels show deep NUV 
imaging with the Swifl/UVOT instrument. The bottom panels show optical/IR images from archival data, as labelled. 



regions to structures generated by th e superposition of clustered, 
unresolved star forming sources (see IPleuss et al.ll200ol for a de- 
tailed analysis of the effects of poorer spatial resolution). The spa- 
tial resolution (2.5 arcsec), pixel size (0.50 arcsec) and galaxy scale 
(74 pc/arcsec) of the images led us to choose the following radii 
for the median kernel: R={5, 10, 15,20,25,30} pixels, which map 
into a physical scale between 0.18 and 1.11 kpc. Th e smallest op- 
tion m aps scales over the largest sizes of HII regions dScoville et al.l 
1200 lb . We combine all those images along with the original one to 
determine, pixel by pixel, the background as the minimum pixel 
value among these images. Fig.[3]compares the UVW2 image with 
the background-subtracted one following this methodology. The se- 
lection criteria applied to SEXtractor, impose a detection threshold 
at a minimum of 9 connected pixels with a count rate 2a above 
the background. The final catalogue comprises 787 sources out to 
R=3h. Fig.|4]shows the radial trend of the number density of NUV- 
detected sources, and those with available Ha photometry - which 
amount to 68% of the total sample. Notice the level in the num- 
ber counts around R=2h (see also inset of Fig.|4j, revealing a shaip 
truncation of the UV light tracing the youngest stellar populations 
in NGC4321. The brightest Ha sources, corresponding to stronger 
sites of star formation, are preferentially located in the central re- 
gions. We will discuss this issue in more detail in §4. The posi- 
tions of the detected sources are used for the photometric analy- 
sis. Note the photometry is performed on the original image, not 
the background-subtracted one. Hence, the photometry can be de- 
scribed as a superposition of a young, star forming region, and an 
older component from the background populations. 



is 91 frames per second. Bright sources will experience a coinci- 
dence loss effect, whereby the true flux is underestimated as all 
photons hitting the detector in a single frame and within a small 
region, ar e counted a s one. This effect has been prope rly charac- 
terized in lPoole etall (2008) and Breeveld et al.l d2010b . and is es- 
pecially important with the optical passbands, where the photon 
count rate is much h igher (as opposed to the UV, see Fig. 1 of 
iHoversten et alj|201 ll . to illustrate this point). 

In our case, optical photometry is available from the Sloan 
Digital Sky Survey (see [J3](, so only the UV data from UVOT 
is used. Fo r the U V images we apply the corrections defined in 
IPoole et al] d2008l) . Our photometric measurements are performed 
within a R=3 arcsec aperture (see §3). For the integration times of 
the final images, the Swift/UVOT exposure time calculate^ gives 
point-source detections at S/N~ 5 of UVW2=23.6, UVM2=22.80 
and UVW1=22.7 (AB), in good agreement with our photometric 
data: Using the detection of a 5a fluctuation above the background 
within the extent of a PSF FWHM, we obtain limiting magnitudes 
in the final images of UVW2=23.5, UVM2=23.2, and UVW1=22.9 
(AB). Fig. [5] translates the UVW2=23.5 limit into a stellar mass 
limit as a function of age - according to the simple stell a r pop- 
ulations from the 2009 versio n of the iBruzual & Charlotl d2003b 
models for a IChabrien d2003l) Initial Mass Function. Both dust- 
less and a dusty, E(B-V)=0.2 mag model are considered. The inset 
shows the distribution of UVW2 sources, with a vertical line mark- 
ing the S/N=5 detection limit. A comparison of photometric mea- 
surements between GALEX/NW (see Fig.© and UVOT/UVM2 - 
which is the filter that tracks best the throughput of the NUV pass- 



Since UVOT is a photon counter, its photometric accuracy for 
brighter sources is limited by the sampling rate of the CCD, which 
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Table 2. Photometry of UVW2 selected sources in NGC4321 1 : SwfWVOT and SDSS (optical). 
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band - gives good agreement, with an RMS scatter for the differ- 
ence NUV— UVM2 around 0.18 mag (taking sources brighter than 
UVM2=22AB). This scatter is mainly dominated by the fact that 
the measurements are taken over a crowded field. 



3 A PANCHROMATIC VIEW OF NGC 4321 

In addition to NUV data, we retr ieved the optical i mages from 
the Sloan Digital Sky Survey/DR8 dAihara et alfeOllh . NGC 4321 
falls at the intersection of four different SDSS fields (run=3631, 
field=456,457; and run=4381, field=99,100). We used the mo- 
saic making facility from the SDSS-III Service Archive Server 
(http://data.sdss3.org/mosaics I to retrieve a properly calibrated im- 
age in all five bands (u, g, r, i, z) and registered the images. In ad- 
dition, imaging is available from Spitzer/IRAC, and KPNO (Ha). 
The reduced, calibrated images are taken fr om the Fifth data deliv - 
ery (April 2007) of the SINGS catalogue^ dKennicutt etaOl 20031) . 
All the optical and IR images are re-sampled to a common size and 
pixel scale, given by the UVW2 frame, which is the image used for 
the selection of targets. We use the IRAF task wregister to ensure all 
bands are aligned properly and with the UVOT original pixel size 
(0.502 arcsec). We compared the photometric estimates before and 
after this transformation, to find no significant systematic change. 
Fig.rjJshows some of these frames, from the UV to the IR. 

The photometry for all images (UV, optical and IR) is per- 
formed with our own code, described below, motivated by the need 
for a consistent background subtraction and error estimate. For each 
position, determined by the SExtractor run, we measure the flux 
within a R=3 arcsec aperture. At every pixel the background is es- 
timated within a circular region with radius R=20 arcsec, where all 
pixels already allocated to a source by SExtractor (i.e. on the seg- 
mentation map) are removed. This distribution is then clipped at 



the 3<t level and the median of the resulting set is taken as the 
background. The IRAC measurements include standard aperture 
corrections as prescribed in the instrument handbook^- We note 
that the standard calibration of the UVOT instrument is done on 
photomet ric measurements of point sources wit hin a R=5 arcsec 
aperture jPoole et"al]|200 it iBreeveld et al.ll2010l) . Given the spa- 
tial resolution of the camera, this choice of aperture is optimal for 
the photometry of Gamma Ray Bursts - the main target of UVOT 
science. However, it does not take advantage of the spatial resolu- 
tion in this case. Hence, we chose a smaller radius for the measure- 
ments (R=3 arcsec), and applied an aperture correction to align with 
the standard calibration. This correction is obtained by performing 
aperture photometry on a number of isolated stars in the frame. 
We selected five bright unresolved targets, for which the photom- 
etry is compared between a R=3 arcsec and a R=5 arcsec aperture. 
This result is extended to the SDSS, and NOAO/Ha frames, to cor- 
rect for the difference in PSF between the different images. For the 
Spifzer/IRAC photometry, we already make use of the published 
aperture corrections, as explained above. The corrections for the 
UVW2, UVM2, UVW1 filters are 0.150, 0.215, and 0.168 mag, 
respectively. In contrast, the aperture corrections for the SDSS im- 
ages are 0.030, 0.045, 0.092, 0.073, and 0.052 mag, in the u, g, 
r, i, and z bands, respectively. Finally, the NOAO Ha image needs 
a correction between these apertures of AlogHa = 0.126 dex. 
Tables 2 and 3 show an excerpt of the aperture corrected pho- 
tometry for the sources in all Swift/UVOT, SDSS, KPNO/Ha and 
Spitzer/IRAC bands, with the error bars quoted at the la level. The 
photometry in the tables is corrected for ( Milky Way) foreg round 
extinction according to the reddening law o f lFitzpatrickld 19991) for 
E(B-V) =0.026 towards the di rection of NGC 4321, taken from the 
maps of ISchlegel et alj ( 1998). 



3 |http : / /sings . stsci . edu| 4 http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/iracinstrumenthandbook/ 
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Figure 3. Illustration of our background subtraction method. The scale is the same for both images. The circles are the apertures used for the selection of the 
background at a given point (see text for details). 



Table 3. Photometry of UVW2 selected sources in NGC4321 1 : Spitzer/IRAC and Ha. 



ID R/h [3.6/nm] [4.5/jm] [5.8/nm] [8/j.m] log Hot 



001 


0.072 


13 


.19±0. 


09 


13.06±0. 


10 


12 


,27±0 


.07 


11 


05±0. 


07 


2 


.31 ±0.04 


002 


0.103 


13 


.02±0. 


12 


12.91±0. 


11 


12 


.05±0 


.14 


10 


81±0. 


15 


2 


06±0.07 


003 


0.104 


13 


,02±0. 


10 


12.89±0. 


09 


12 


.23±0. 


.07 


11 


02±0. 


05 


2 


,44±0.03 


004 


0.116 


12 


,68±0. 


07 


12.59±0. 


08 


11 


88±0 


.06 


10 


68±0. 


06 


2 


.16±0.05 


005 


0.212 


15 


.50±0. 


12 


15.47±0. 


10 


15 


,08±0 


.09 


13 


99±0. 


07 






006 


0.254 


16 


.96±0. 


13 


16.97±0. 


10 


16 


.95±0. 


.12 


16 


60±0. 


22 





26±0.08 


007 


0.276 


15 


.79±0. 


08 


15.60±0. 


06 


14 


.76±0. 


.09 


13 


49±0. 


06 


1 


.32±0.07 


008 


0.309 


17 


.25±0. 


07 


17.30±0. 


08 


16 


.90±0 


.07 


15 


81±0. 


09 





.24±0.33 


009 


0.312 


16 


.34±0. 


09 


16.29±0. 


07 


15 


.40±0 


.07 


14 


23±0. 


09 





,94±0.07 


010 


0.320 


16 


.58±0. 


14 


16.44±0. 


13 


15 


.98±0 


.10 


14 


86±0. 


11 





,86±0.08 


011 


0.327 


16 


.46±0. 


10 


16.33±0. 


08 


15 


.12±0. 


.08 


14 


00±0. 


08 


1 


.25±0.06 


012 


0.330 


17 


.38±0. 


09 


17.44±0. 


08 


17 


,11±0 


.10 


15 


90±0. 


09 





.34±0.07 


013 


0.331 


16 


80±0. 


08 


16.75±0. 


09 


16 


.22±0. 


.08 


15 


08±0. 


07 





.87±0.03 


014 


0.349 


17 


.03±0. 


13 


16.79±0. 


11 


15 


.61 ±0 


.12 


14 


.31±0. 


09 


1 


16±0.03 


015 


0.352 


17 


.80±0. 


14 


17.85±0. 


09 


17 


.42±0. 


.11 


16 


24±0. 


11 


-0 


.04±0.28 


016 


0.358 


17 


.82±0. 


13 


17.80±0. 


14 


17 


.19±0. 


.15 


15 


76±0. 


17 


-0 


.00±0.22 


017 


0.361 


17 


.42±0. 


14 


17.42±0. 


16 


18 


.18±0. 


.88 


20 


61±1. 


19 


-0 


.96±0.08 


018 


0.370 


17 


.97±0. 


11 


17.93±0. 


15 


18 


.67±0. 


.62 











.70±0.06 


019 


0.387 


18 


.35±0. 


23 


18.40±0. 


17 


19 


.34±0. 


.68 











.74±0.04 


020 


0.391 


17 


.81±0. 


15 


17.46±0. 


16 


16 


.77±0. 


.16 


15 


.77±0. 


59 





.85±0.17 


021 


0.406 


17 


,33±0. 


10 


17.03±0. 


06 


15 


.40±0 


.09 


14 


10±0. 


08 


1 


.25±0.08 


022 


0.415 












18 


.78±0. 


.43 


16 


55±0. 


45 






023 


0.417 


18 


.56±0. 


22 


18.81±0. 


28 


















024 


0.428 


17 


.75±0. 


12 


17.56±0. 


12 


16 


.45±0. 


.14 


15 


22±0. 


17 





.82±0.06 



1 Photometry for the IRAC bands in the AB system, uncertainties quoted at the 1<t level (random errors 

only). This is an excerpt of the full table that can be accessed in the electronic version of this paper. 
T measured as log 10 [L(Ha)/10 37 erg/s]. 



4 MODELLING THE STELLAR POPULATIONS 

Fig. [6] shows the NUV and optical colours of the UVW2-selected 
sources. A number of colour-colour diagrams are shown, involving 
UVOT and SDSS photometry. The observations (dots) match the 
predictions for a set of simple stellar populations (SSP). All mod- 
els have solar metallicity. The thick soli d orange lines correspond t o 
dustles s SSPs from the latest models of lBruzual & Charloj d2003l) . 
with a lChabrierl d2003h initial mass function (IMF). To illustrate 



the effect of using a different population synthesis model, or a dif- 
ferent IMF, we also show in the figure the colour tracks from the 
Starbu rst99 dustless models (Leitherer et al. 2010}) with a|Salpeter 
( 1955) IMF (blue dashed lines). For reference, crosses mark stellar 
ages of 1, 100, 500 and 900 Myr in the latter. We also include, as a 
thin solid red line, the effect of a E(B-V)=0.2 mag dust screen on 
the Braz ual & Chariot mo dels, according to a Milky Way redden- 
ing law dFitzpatrick|[l99"3) . The use of this reddening law is justi- 
fied by the type of dust present in NGC4321 (see e.g. lDraine et al.l 
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Figure 4. Radial profile of the surface density of UVW2 sources (thick 
black line), including the density of counterparts in the Ha images (red). 
The surface density of bright Ha sources (L(Ha)> 10 37 erg/s) is shown 
as a dashed red line. The error bar give the Poisson fluctuations. The inset 
shows the extent of the UVW2 disc, with projected circles - according to 
the orientation of the disc - representing the R=h, 2h, 3h radii, where h is 
the characteristic scale-length measured in the V-band jKodaiiaetalJl98d) 

120071) . The model predictions clearly do not favour large amounts 
of dust, as expected from a UV-based source detection. Notice that 
on the colour-colour diagrams the dust content of the sources is 
found within a small range of E(B-V) ,$ 0.2 mag. We will quan- 
tify this statement in some detail below. Fig. [6] also reveals that 
the age of these sources is never older than about 1 Gyr, although 
we have to take care of the presence of background populations of 
older stars. Note the UV colours are in better agreement with the 
SSP models, whereas the optical colours are more dispersed. This 
is caused by the contribution at redder wavelengths from older pop- 
ulations within the apertures. We explore this issue in more detail 
in the next figure. 

Fig. UJ shows the trend of the Ha luminosities with respect to 
UVW2 magnitude, for those sources with a detection in the Ha 
image (68% of the total sample). Ha emission in star forming re- 
gions originates from gas excited by ionising photons from OB 
stars, i.e. more massive than ~ 15M0. In contrast, NUV emission 
will also originate from the photospheres of lower mass stars, down 
to ~3M0 . Hence, the ratio between Ha and NUV light works as a 
stellar clock over the first few million years of the nascent popula- 
tion (assuming a universal initial mass function, a s expected for a 
bright galaxy, see iHoversten & Glazebrookll2008l) . The solid lines 
in the top panel of Fig.[7JcoiTespond to the model predictions from 
Starburst99 for a simp le stellar population at solar metallicity, with 
a Kroupa IMF (see lLeitherer et ai]|201(l and references therein). 
The Kroupa IMF is very similar to the Chabrier IMF used below 
for the analysis of the photometry: differences between them at the 
level of the data used are indistinguishable and do not affect our 
conclusions. The three solid lines represent three different stellar 
masses, from 10 3 to 10 5 Mq. As age progresses, the Ha flux de- 
creases: the crosses, from top to bottom, mark the ages of 1, 5 and 



Figure 5. The stellar mass is shown with respect to the age of a simple 
stellar population for the limiting magnitude of UVW2=23.5 (at S/N~5). 
Dustless (solid) and a dusty case (dashed) is considered, as labelled. The 
inset shows the histogram of UVW2 detections, with the vertical line locat- 
ing the S/N~ 5 case. 

10 Myr. The star forming sources are coded according to radial 
position. To highlight the differences with respect to galactocen- 
tric distance, we show those sources with R<h (solid red dots) or 
R>2h (open blue dots). The radial trend is evident, with the central 
sources being more m assive, a result consistent with the findings of 
iGoddard et all JioT ]]). Notice the orientation of the "dust vector" 
for E(B-V)=0.2 mag. If dust were an important factor in the inner 
regions (i.e. the solid red dots), the segregation between central and 
outer sources in this graph would be more prominent, as a redden- 
ing correction would shift the points corresponding to the central 
regions towards the right of the figure, therefore increasing the sep- 
aration between inner and outer sources. The ages of all sources 
are mostly younger than 10 Myr, as expected for an Ha detection. 
The bottom panels show the optical colours, with the solid lines 
representing the Starburst99 predictions for the same models as in 
the top panels. Differences with respect to the radial distance are 
less evident, with a tail towards bluer colours in the outer sources. 
Those sources with colours redder than those from the simple stel- 
lar populations reflect the contribution from an underlying older 
component, which will dominate optical and NIR light. This figure 
illustrates the shortcomings of simple stellar populations in describ- 
ing the photometric data. 

The simplest model to overcome this problem involves a su- 
perposition of two SSPs. These models have been applied for the 
detection of residual star form ation over the old populations found 
in ear ly-type galaxies (see e.g. lFerreras & S ilk 2000; Kavir aTet al.l 
120071) . The models we use in this paper are described by five free 
parameters comprising the age of the old and young components 
(to and ty, respectively), the metallicity of both populations (as- 
sumed to be the same), the mass fraction in young stars (fy) and 
the dust reddening affecting the young component (Ey), whereas 
the old component is considered to be dustless. In order to deter- 
mine robust best fits, we opted for a full search of this five dimen- 



© 2012 RAS, MNRAS OOO.mmi 



8 Ferreras et al. 



10 20 30 40 



1 

0.5 

I ° 
-0.5 

-1 
2 

m 1 



-1 
1 

S 0.5 
% 
5 -0.5 



-_ 1 11 1 _ 


-_ 1 1 1 1 1 1 1 1 1 1 1^ 




1 1 1 II 1 

\jf : 

~ 1 1 1 II 1 " 


1— 1 1 1 1 1 1— 1 ! 1 1 1-^ 

q — i — h- — i — h- 1 — |— i — i — i — | — f 


1 : , ,' 1 ■ , , , , 1 , , , , 1 , , : 


: , ,1 , ' , ■, 1 , , , 1 , : 



_ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 l_ 

p ^ oR>2h 


r j — 


:' O , 

<-Tt^^^telLv«# • * "*. 

;°<> * ^ ; 
i , , , r» ,• , , i , , , i , , , r 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

T~l 

l_L . 
1 1 

n 
u 

ll.nl..>. I....I." 



0.5 1 

UVW2-UVW1 



2 4 

UVW2-g 



24 



22 20 IE 

UVW2 (AB) 



16 10 20 30 40 



Figure 6. Colour-colour diagrams of the star forming regions in NGC 4321. 
The observed magnitudes have already been corrected for foreground ex- 
tinction (E(B-V)=0.026). The models correspond to dustless CB09 (thick 
solid line, orange), SB99 (dashed line, blue) and CB09 with a E(B- 
V)=0.2 mag dust screen (thin solid, red). See text for details. The SB99 
models in the bottom panels are marked by crosses for ages (from left to 
right) of 1, 100, 500 and 900 Myr. 



sional parameter space. The grid of parameters is shown in Table[4] 
involving over 100 million models. We compute the flux from the 
SSPs, integrated over the filter response functions, and store them 
in a lookup table, to maximise the computing speed of the fitting al- 
gorithm, which is also parallelized to reduce wall-clock computing 
time. The effect of dust on the magnitudes is modelled by a fourth 
order polynomial in E(B-V), where the coefficients are computed 
from a grid of 16,384 SSP models over a range of ages, metal- 
licities and E(B-V), using a lFitzpatrickl ( ll999h reddening law. For 
each choice of parameters, the photometric data from the UVOT, 
SDSS and IRAC channels 1 (central wavelength A c = 3.6/im) 
and 2 (A c = 4.5/xm) are used to define a likelihood. Following 
a Bayesian approach, we apply a Gaussian prior on the metallicity, 
with mean log Z/Zq — —0.1 and standard deviation 0.3 dex. The 
choice of mean is based on the observat ions of the gas-phase metal- 
licity of Ml 00 jMoustakas et alJ f201QD. However, we note that by 
removing the prior altogether, the age distribution of the old and 
young components is not affected within the scatter found among 
the sources. 

Fig.[8]shows the distribution of the ages - total mass weighted 
average (top) or age of the youngest component (bottom) - and dust 
reddening for all sources with respect to UV-optical colour or UV 
flux. The big dots track the median of the distribution within bins 
keeping equal number of sources per bin, with the error bars rep- 
resenting the RMS scatter. Solid (open) dots correspond to sources 
at radial distance from the centre R<h (R>2h). As expected, UV- 
optical colour tracks the age of the young component and no sig- 
nificant trend is found with respect to the dust, obtaining a value 
for the reddening around E(B-V) 0.2 mag. The distribution of 
metallicities has a mean of log Z/Zq — —0.18 - as expected, 
slightly lower than the measured gas-phase metallicity - and an 



Figure 7. The Ha line luminosities and optical colours of the star form- 
ing regions are shown with respect to UVW2 magnitudes. The star form- 
ing regions are split with respect to galactocent ric distance, as labelled. 
The lines track the predictions from Starburst99 (Leitherer et al. 2010) of 
a simple stellar population with solar metallicity, and stellar masses be- 
tween 10 3 and 10 5 Mq. The crosses mark (from top to bottom) ages of 
1, 5 and 10 Myr. The panels on the right show the distribution of inner 
(dashed) and outer (solid) star forming regions. The arrows in both panels 
correspond to a dust attenuation of E(B-V)=0.2 mag using a Milky Way 
reddening law. The dashed line in the top panel marks a limiting line lu- 
minosity logL^ (erg/s)= 37.9 — 0.4(UVW2ab — 20), used as a proxy to 
segregate Ha detected sources with respect to age (see Fig. II U . 



RMS scatter of 0.08 dex. There is no significant radial trend in 
metallicity, although the RMS scatter is higher at R>2h (0. 12 dex). 
Note that the ages of the young component found here are signifi- 
cantly older than those estimated from the Ha luminosities (Fig. [7}, 
reflecting the co mplex mixture of populations within a resolution 
element (see e.g. lPleuss et al.ll2000h . Within a UVOT detection, we 
expect to have a wide range of stellar populations, including dust 
obscured regions, Ha emitting sources, and more mature UV-bright 
areas. An analysis based on the broadband photometry will be in- 
herently biased towards the last component. The presence of dust 
obscured regions can be illustrated by the stacking of sources into 
two representative UV-to-IR spectral energy distributions (SEDs). 
Fig.|9]shows the photometric SEDs of the youngest sources (using 
the UVW2-g colour as a proxy for age) have an IR excess in the 
IRAC bands, which reflects the contribution from dust-obscured re- 
gions that are possibly too faint in the UV-optical range to have an 
effect on the analysis of the UVOT/SDSS photometry. However, 
this figure gives further support to the use of UVW2-g as a proxy 
for the stellar ages of these sources. Notice in Fig. [8] the lack of 
any trend in the properties of the stellar populations with respect to 
radial distance, except for the fact that the outer sources are fainter 
in UV (see also Fig. [7}, confirming that the radial change refers to 
the stellar mass content of the sources. The young component has 
a distribution of ages in the range 10-100 Myr. We will use this 
information in the next section to explore in detail the distribution 
of UV sources with respect to the spiral arms. 
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Figure 8. Distribution of ages and dust reddening of the UV sources ac- 
cording to the two-burst model (see text for details). The solid (open) dots 
show median values for subsamples at R<h (R>2h), with the binning cho- 
sen at equal number of sources per bin. The error bars correspond to the 
RMS scatter of individual sources within each bin. The top panels show the 
mass-weighted average age, whereas the bottom panels correspond to the 
age of the youngest component in the two-burst model adopted here (see 
text for details). 



Figure 9. Spectral Energy Distribution of a subsample, selected according 
to UVW2-g colour as labelled. The lines and shaded areas map the average 
and standard deviation of all sources within the colour selection, shown as 
the flux through an "X" band representing the filters from UVOT to IRAC, 
as shown in the horizontal axis, measured as a ratio with respect to the 
SDSS r flux. The blue sources (top panel) have a clear excess towards longer 
wavelengths, possibly from PAH dust emission. 



5 PUTTING DENSITY WAVE THEORY TO THE TEST 

The properties of the underlying stellar populations of the star 
forming sources provide valuable constraints on the origin of 
spiral arms. In the previous section, we showed that our NUV- 
detected sources span a wide range of stellar ages (a few hun- 
dreds of Myr) and are distributed over a large range of galacto- 
centric radii. Therefor e, they represent an ideal tracer for the "off- 
set method" (see e.g. lEgusa et ai] l2009l) . to look for variations in 
the distribution of different star formation tracers around the spi- 
ral arms. Standard density-wave theory predicts a systematic trend 
in the distrib ution of stel l ar age s inside and outside of the coro- 
tation radius. iFovle et al.l fcOl lh used UV images from GALEX 
and IRAC/3.6/im emission to trace older populations. However, 
NGC 432 1 is not in their sample and the spatial resolution of UVOT 
is superior to GALEX. 

In order to analyse the properties of the star forming sources 
with respect to the spiral arms, we fit three portions of the spiral 
structure of NGC 4321 in regions where the arms can be traced. 
Fig.QJJshows the selection of the UV sources along the spiral arms. 
We use the UVW2 and the Ha images to manually select a number 
of points along the spiral arms. The solid lines tracing the arms in 
Fig. [10] are obtained from a linear interpolation of these sources 
in polar coordinates. The crosses, triangles and circles locate the 
UVW2 sources that reside within ±2 kpc of the spiral arms. 

Fig. [TJJ (left) shows a sketch of the changes expected in the 
distribution of the spiral arm offsets of the star forming sources 
with respect to age and galactocentric distance. The solid blue and 
dashed red distributions illustrate the trend of the offsets between 
a very young component (e.g. an HII region) and an older one 



(e.g. an UV star forming source), respectively. The telltale sign 
of a density wave with constant pattern speed would be the shift 
of the older component with respect to the young one inside and 
outside corotation. We emphasize here that we only focus on rela- 
tive offsets between these two components, noting that there would 
be an additional, smaller, offset between the HI I regio ns and the 
molecular gas, as found in, e.g. Tamburro et all (l2008h . The his- 
tograms on the right hand side of Fig. II II show observed distri- 
butions of various properties of the sources with respect to the 
distance to the spiral arm for two choices of ra dial separation, at 
either side of the corotation radius suggested by ICanzian & Alien! 
( 1997) (top panels: R>1.5h; bottom panels: R<h). They compared 
the two dimensional observed velocity field with what the first- 
order linear spiral density wave theory predicts, and measured a 
corotation radius R cor = 98 ± lOarcsec (shown for reference 
in Fig. 110b. We note that this estimate, located around R~ 1.5h 
falls in the middle of our distribution of sources. We will adopt 
this value for reference in our analysis, although we note that 
other studies give higher results: Sempere et all (1995) propose a 
corotation radius in the range 8-11 kpc, corresponding t o 100 - 
150 arcsec. The pattern speed measured bv lRand & Wallinl d2004l) , 
using the Tremai ne- Weinberg me thod, imp lies a corotatio n radius 
of ~ 2arcmin fegusa et alj|2009t) , while lOev et al.l J2003I) suggest 
a corotation radius of 2.6 arcmin. However, these methods may not 
be applicable if the spi ral arm pattern spe ed changes as a function 
of radius. For instance. iMeidt et all d2009l) applied their developed 
version of the Tremaine- Weinberg method and demonstrated that 
some spiral galaxies show a radial variation of the spiral pattern 
speed. Unfortunately, NGC 4321 is not in their sample. 

The left and central panels are model-independent compar- 



© 2012 RAS, MNRAS OOO.mfTH 



10 Ferreras et al. 



T- 1 | 1 I— — 1 | 1 I 1 | . I .1 

6 ... . ' . 




J I I !_, I -L 

2 

ARA (arcmin) 



Figure 10. The Ha sources are used to define three portions of the spiral 
arm structure where we can study the distribution of star forming sources. 
The contour map shows, as reference, the UVW2 image, which tracks a 
wider range of stellar ages. The dashed projected circles have radii h, 3h/2, 
and 2h, where h is the disc s cale-length. The thick so lid circle locates the 
corotation radius according to lCanzian & Allenl jl997h . The crosses, trian- 
gles and circles are the UVW2 sources located within ±2 kpc of the defini- 
tion of the arms. 



isons of NUV-detected sources with respect to their Ha luminos- 
ity or UV-optical colour, respectively. The Ha selection is split 
at a reference value for the luminosity of the line denned by 
logL°(erg/s)= 37.9 - 0.4(UVW2 AB - 20), motivated by the lo- 
cation of the Ha sources in the top panel of Fig. [7] Hence L„ acts 
as a Ha-based tracer of stellar age. The middle panel is motivated 
by the broadband photometric analysis (see Figs.[8]and|9) but does 
not rely on any modelling. Finally, the rightmost panels separate the 
sample according to the ages of the young components according to 
this analysis, and described in the previous section. No significant 
offsets are found in any of these cases. 

The lack of an offset is clearly inconsistent with a solidly 
rotating spiral arm pa ttern as predicted by density wave theory. 
iDobbs & Pringlel ( l2010l) show a predicted distribution of young star 
clusters with different ages from their numerical simulations with a 
fixed spiral arm potential rotating solidly and mimicking a long- 
lived spiral arm. Their fig. 4 provides a quantitative plot of the 
sketches we outlined on the left of our Fig. QT] and they consider 
alternative mechanisms for the formation of spiral arms that would 
give rise to similar distribut ions as our observed histograms. Inter- 
estingly, [Egusaet al. (2009) also found no offset between CO and 
Ha spiral arms in NGC4321, although we note that their analy- 
sis covered a narrower radial range (R< 1.4', roughly one scale 
length). In combination with our analysis, we can therefore state 
that in NGC4321 no offset is found out to 2 scale lengths over a 
wide range of ages, from molecular clouds to about 100 Myr young 
stars. This is a stri king result and agrees with what is also found in 
iFovle et al.1 bOllh for other galaxies. Our result provides further 
evidence against classic spiral density wave theory. 

The distributions found in Fig.QT]can be naturally explained 



if the spiral arm is "corotating" with the stars and the gas, i.e. if 
the spiral pattern speed decreases with galactocentri c radius. Re- 
cent numerical simulation studies dWada et all 1201 ll : iGrand et"ai1 
120121) claim that such corotating spiral arms are more naturally de- 
veloped in numerical simulations. Indeed, it is rather difficult to 
reproduce the spi ral density wave theory in numerical simulations 
with self-gravity (Sellwood 201 1). Such corotating spiral arms suf- 
fer from the winding problem, but are found to be disrupted within 
a time-scale of a few hundreds of M yr (see e.g. iFuiii et al]|201 ll : 
lOuillen et alj|201 H ; Isellwood 2011, and references therein) and a 
snapshot of the evolving g alaxy at an y time always shows some 
spiral arm structure. Wada et all (1201 lb also demonstrate that in 
such corotating spiral arms gas falls into the spiral arms both from 
behind and the front of the spiral arm and collides, inducing star 
formation in the arm. This scenario is consistent with the observed 
distribution of star forming sources presented here. 



6 CONCLUSIONS 

The resolution of the Swiff/UVOT camera makes it an optimal in- 
strument to explore the young star forming sources in nearby late- 
type galaxies. In combination with optical photometry, it is pos- 
sible to disentangle the contribution from the more diffuse, older 
component, giving an estimate of the age and dust content of the 
young stars responsible for the NUV emission. In this study, we 
combine deep UVOT imaging in the NUV with optical and IR 
archival imaging to study the properties of 787 star forming re- 
gions in the spiral structure of NGC 4321, out to R~3h. We find a 
strong tr end of brighter Ha sources in the central regions (Fig. [7] 
see also Goddard et alj|201l[) . whereas the ages, or the dust con- 
tent of the star forming regions do not present a measurable cor- 
relation with galactocentric distance (Fig. [HJ, suggesting that the 
strong trend in Ha flux is caused by a larger number of massive 
star forming sources in the central regions (assuming a universal 
stellar initial mass function). The presence of dust is revealed when 
stacking up the photometric SEDs of the sources according to their 
NUV-optical colour, with significant excess at 8/j.m for the bluest 
(and youngest) ones (Fig. |9). However, a quantitative analysis of 
the dust reddening from the photometric data gives rather low "spa- 
tial resolution averaged" values of E(B- V) ,$ 0.2 mag. The main 
outcome of the comparison with stellar populations is that one can 
use UVW2-g colours as a proxy of the age of the younger compo- 
nent, ranging from a few Myr to a few hundred Myr (see Fig. |8}. 

The properties of these sources enable us to study the effect 
of the spiral structure, a powerful test of pattern formation in disc 
galaxies. The standard scenario of density wave theory with a con- 
stant pattern speed results in an offset with respect to age for the 
distribution of distances to the spiral arms as one moves from the 
central regions - where the gas moves faster than the pattern - to 
the outer disc (see Fig. [TJJ left). The advantage of UV imaging 
over studies based on comparisons between the gas component and 
HII/24/im regions is the wider time baseline spanned by UV light, 
roughly of order 100 Myr. No significant differences are found in 
the distribution of these sources, giving further negative evidence 
for density wave spirals. Our fin dings reject the possibility that the 
results from Eg usa et aTl (|2009) - giving no offsets between CO 
and Ha over a narrower range (R< 1.4') - can be explained by a 
radial coincidence with corotation. Instead, a scenario with short- 
lived corotating spiral arms is compatible with the data. Such struc- 
tures have been found in numerical simulations of spiral galaxies 
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Figure 11. Left: Sketch of the expected shifts in the distribution of stellar populations according to density wave theory. A spiral pattern with constant angular 
speed - represented by the vertical line - interacts in different ways with the rotating gas depending on whether the gas is located inside or outside the 
corotation radius Rcor- On the left panels, a comparison of very recent star formation - as tracked by Ha (solid blue curves) - with older stars - from UV 
light (dashed red curves) betrays the presence of a pattern in solid body rotation. Right: Comparison between the distribution of UVW2-selected sources split 
according to Ha line luminosity (L(Ha), leftmost panels), UV-optical colour (UVW2-g, middle panels) or the age of the young component (ty, rightmost 
panels), estimated from the modelling of the NUV/Optical/IR photometry. The threshold in Ha luminosity, , is defined with respect to the UVW2 flux 
m Fig.|7] Only sources within ±2 kpc from the three spiral arms defined in figure ITol are included in the histograms, split according to radial distance to the 
centre, measured in disc scale lengths. 



Table 4. Grid of 2-population models used in the analysis 



Parameter 


Young 


SSP Old SSP 


Models 


Age (Gyr) 


-3 < log 


t Y < < logto < +1.12 


64 X 64 


Metallicity (Z Q ) 




-2 < logZ < +0.3 


24 


Dust Reddening 


< E B _ 


_v < 1 dustless 


32 


Mass fraction 


0</y 


< 0.8 1 - / Y 


32 



TOTAL 1.007 x 10 8 



(see e.g. iGrand et al.ll2012l) . implying a fundamental change of our 
views from the traditional concept of long-lived spiral arms. 
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